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We study the Ising model in a hierarchical small-world network by renormalization group analysis, 
and find a phase transition between an ordered phase and a critical phase, which is driven by the 
coupling strength of the shortcut edges. Unlike ordinary phase transitions, which are related to 
unstable renormalization fixed points (FPs), the singularity in the ordered phase of the present 
model is governed by the FP that coincides with the stable FP of the ordered phase. The weak 
stability of the FP yields peculiar criticalities including logarithmic behavior. On the other hand, 
the critical phase is related to a nontrivial FP, which depends on the coupling strength and is 
continuously connected to the ordered FP at the transition point. We show that this continuity 
indicates the existence of a finite correlation-length- like quantity inside the critical phase, which 
diverges upon approaching the transition point. 

PACS numbers: 64.60.aq,75.10.Hk,64.60.ae,89.75.Da 



Recently, various physical phenomena in non- 
Euclidean graphs have been studied especially in the con- 
text of complex networks, and their properties have beer 
found to be beyond the scope of the conventional theorj 
for Euclidean graphs 1 . Of particular interest are sys- 
tems regarded as infinite dimensional in the sense thai 
the equidistant surface S r of radius r grows exponen- 
tially as S r oc e VdT with a positive constant yd, which i; 
faster than any power function r d_1 as in d-dimcnsiona 
Euclidean graphs. Typical examples are trees and hy- 
perbolic lattices 2 . Remarkably, such infinite-dimensiona 
systems often exhibit the critical phases in which tht 
(nonlinear) susceptibility diverges 3-5 . Although the crit- 
ical phase is also observed in some Euclidean systems 
e.g., the quasi- long-range ordered phase in the two- 
dimensional XY model 6,7 , the critical phases in infinite 
dimensional systems are considered to be due to rathci 
geometrical effects. Indeed, the exponential growth or 
a graph admits the divergence of the susceptibility x, 
which is calculated by the integral of the two-point cor- 
relation function 20 , even with finite correlation length £ 8 ; 
for | > l/y d 




X oc 



In N/y d 



drS r e- r/ t oc N*, tp=l- l/y4 (1) 



diverges in the thermodynamic limit N — > oo. Here we 
set the upper bound In N/yd to ensure x oc iV for ^ = oo. 
A critical phase, if it exists, lies between a disordered 
phase with £ < l/yd and an ordered phase with £ = oo. 
Such a phase with a fractal exponent < ip < 1 is actu- 
ally observed in the percolation transitions in enhanced 



FIG. I: (a) Hierarchical small-world network with eight 
nodes. The periodic boundary condition is imposed on the 
horizontal direction, (b) One-step decimation of the central 
spin in the partial sum of partition function. This is an inverse 
procedure to grow the graph in (a). 



trees 5 , hyperbolic lattices 9 , hierarchical graphs 10,11 , and 
growing random networks 12 . 

The property of the phase transitions between a crit- 
ical phase and an ordered phase is an interesting issue. 
Quite recently, some models in the simple hierarchical 
network shown in Fig. 1(a) were investigated to exam- 
ine these transitions. This network is very useful be- 
cause rigorous real-space renormalization is possible for 
various models in the simplest way. Furthermore vari- 
ous types of phase transition are observed depending on 
the model used, e.g., a discontinuous transition of the 
bond percolation model 11 , equivalent to the 1-state Potts 
model 13,14 , and continuous transitions with a power-law 
singularity (PLS) or an essential singularity (ES) for the 
q-state Potts model with q > 3 15 . These are observed in 
other graphs 3-5,16,17 . Thus this hierarchical network is 
a good stage to investigate what determines the type of 
phase transitions in a systematic way. In particular the 
2-state Potts model, equivalent to the Ising model, stands 
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between q = 1 and 9 = 3, and has special importance to 
understand how the transition class changes. 

In this Rapid Communication, we study the phase 
transition of the 2-state Potts model in the network men- 
tioned above by renormalization group (RG) analysis, 
which reveals a new class of phase transition governed 
by the stable fixed point. Furthermore the singularity 
for q = 2 is found to be special due to the marginal bi- 
furcation of the RG fixed point (FP) between pitchfork 
type for q < 2 and saddle- node type for q > 2. 

We consider the 2-state Potts model in the hierarchical 
small- world network shown in Fig. 1(a). This consists 
of one-dimensional backbone edges (BBEs) and nested 
shortcut edges (SCEs). The spin variable a, taking a 
value or 1, is put on every node of the network. The 
dimensionless energy function is written as 
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where K and J are the coupling constants on BBEs 
and SSEs, respectively, and H{> 0) denotes the exter- 
nal magnetic field coupled with the state 0. 

Let us consider the partial sum of the partition func- 
tion Z — Yii er E l kBT over the states of spins in the 
youngest generation having two neighbors [see Fig. 1 (b)] . 
This is rigorously equivalent to replacing the parameters 
as 

O „K n+ i5 r , ia2 +(H n+1 /2)(5 ai o+S a2 o) 

with the recursion relations: 

g n+1 = g n + 2-(" +1 > [K n + J- + A™ / 2 + A^] , 

K n+1 = - (K n + J) + A^/2 + Al? /2 - A& , (4) 
H n+1 = H n + AW -A%\ (5) 

1(2) 



where g„ 



lnC„, e A ^ 
^ _ 



2(K n + J)+H n 



1. e 



e 2(K n +J) + e H, % and = £ H n + L 

For H = 0, the recursion relation is rewritten as 



Wi+l 



K [fkl (2j 1)] 



(6) 



where we put k n 
by k n+ i = k„ 

k* = 



-K„ 



and j 



-J 



The FP is given 



k* , which is solved as 



0,±y/2t/(l + t), t = 2j-l. 



(7) 



There is a pitchfork bifurcation point (PF BP) at (j, k) = 
(1/2,0); the stable FP is k* = for j < 1/2 and 
k* = y/2t/(l + t) for j > 1/2 (see Fig. 2). The re- 
gion where the flow goes to the nontrivial fixed line is 
regarded as a critical phase because each FP on the line 
represents a self-similar structure like an ordinary criti- 
cal point. It is noteworthy that in this renormalization 




FIG. 2: (color online) Phase diagram in the j vs k plain. The 
shaded region indicates the critical phase. The (blue) solid 
line and (red) dashed line denote stable and unstable fixed 
lines, respectively. The arrows indicate the RG flow of k. 
The inset shows the fixed lines for q=l, 2, 3, 5, and 8. 



procedure the couplings on the SCEs remain the same 
while the ones on the BBEs are updated. The SCE cou- 
plings < Jc 1 prevent K n from converging to zero, 
and therefore the critical phase appears instead of a dis- 
ordered phase. A similar mechanism is observed in the 
RG analysis of the Ising models exhibiting the critical 
phases, e.g., in graphs: the decorated (2,2)-flower 3 and 
the Hanoi network with average degree 5 (HN5) 18 . 

When j increases with k fixed, a transition from the 
ordered phase to the critical phase occurs at j = j c = 
1/2, irrespectively of the value of k. Remarkably, the FP 
at j c corresponds to K = oo(k = 0) as well as in the 
ordered phase j < j c and is stable, so that the critical 
behavior is quite unlike the conventional one governed by 
unstable FP (US FP). This yields a curious singularity, 
as shown below. 

Equation (6) is approximated by a differential equa- 
tion: 



2s r , 



dn 



This is solved as e 2n l*l 



\t\) 



2 „-2 



<3 



(8) 



so \ ■ 



+ 1*1 
n + |t| 



1-1*1 



leading to 



2n\t\ 



1 



l/2n for n\t\ < 1 
|t| e -2n|t| for n | t j > 1 



(9) 



for \t\ <C So- Since s n goes to zero for n — > oo, it is (and 
thus K is) an irrelevant parameter. 

In the lowest order of H n , Eq. (9) is still valid as a so- 
lution of Eq. (4). On the other hand, Eq. (5) is rewritten 
as 



+i 



H n + In ■ 



s 2(K n +J)+H„ 



1 



,2(K n +J) _|_ g 



2H n 
1 + s„ 



O(H 2 n \10) 



Note that we leave s n in the first term; the cross term 
of s n and H n remains and thus this is beyond the linear 
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stability analysis. By putting H n = 2 n e Xn , we obtain 
\t\d£ 1, l-e" 2n l*l 



n-l 

t=n 



"(J 



1 1 

e u\t\ _ { ~ 2 i _ e -2"o|t|' 



where we substitute Eq. (9) into the summation and ap- 
proximate it by the integral. Here we give the lower 
bound of summation no to avoid the divergence of sg 
at £ = 0. Finally we obtain 



H n = He nVd 



Un V 1 - e 2-1*1 \ Itl 1 / 2 for n|t| » 1 ' 

where we use yd = In 2. Note that the factor u n does not 
appear in the linear-order treatment for s n and H n . 

The RG evolution Eq. (12) leads to a scale-invariant 
formula of the free energy as 



g(H,N- 1 ) = e- n y«g(Hu r , 



l ,N~ 



') (13) 



for n > 1, Here t is not included explicitly in the ar- 
guments of g because it is (and thus J is) not a scaling 
field but just an external parameter. The irrelevant scal- 
ing field s n is also omitted. However the effect of s n is 
indirectly included in H n as in Eq. (10) and gives the t 
dependence as in Eq. (12). The first and second order 
derivatives of g with H arc written as 



gniH,^ 1 ) = u n g H (Hu n e nyd ,N~ 



- 1 e nyd ), (14) 



respectively. In the following, we show the behaviors of 
the magnetization to = g^ and the susceptibility x = 
gn 2 in three asymptotic regimes: (i) t < 0, H = and 
N = oo, (ii) t = 0, H > and N = oo, and (iii) t = 0, 
H = and TV < oo. 

(i) For t < 0, H — and N = e nyd oo, we obtain 



to = \t\V 2 g H (0, 1) oc |t| x / 2 , 
X = N\t\ gH *(0,l)<xN\t\, 



(16) 
(17) 



with u n — > It] 1 / 2 . The magnetization shows the PLS as 
if it were an ordinary second order transition. Within 
the linear analysis u n = 1, m does not vanish at t = 0, 
and the transition would look discontinuous. Surpris- 
ingly, x diverges in the whole ordered phase, meaning 
the coexistence of a divergent fluctuation and a nonzero 
order parameter. Such coexistence may be realized by 
spatial segregation into the fluctuating region and the 
ordered region around the root node, the left most node 
in Fig. 1(a). 

(ii) For t -> and N = e nyd , u n -> iiT 1 ! 2 , Eqs. (14) 
and (15) become 



gniH^N- 1 ) 
gn^H^N- 1 ) 



1/2 g H (Hn~^ 2 e nyd ,l), (18) 
1 e nyd g H 2(Hn- 1 / 2 e nyd ,l). (19) 




|t| In H" 



FIG. 3: (color online) Scaling plot of the magnetization in the 
ordered phase. Data with m < 0.20 are used. 



We set n so as to satisfy Hn~ 1 / 2 e nVd = 1, which is ap- 
proximately solved as e nyd ps \n(H~ 1 )/ydH. By substi- 
tuting this into Eqs. (18) and (19), we obtain 



!J,I 



^g H (l,l)<x [HH- 1 )] 



-1/2 



ln(i? 



; Qh 2 1,1 OC — . 

In if" 1 y d H JH V ' ; if 



(20) 



(21) 



(iii) For t = and if = 0, Eqs. (18) and (19) become 



y^ 5ff (0,l)cx(ln7V)- 1/2 , (22) 

(23) 



N N 
9h 2 (0, 1) oc 



In TV 



In TV' 



On the crossover between the two limits, Eqs. (16) and 
(20), we expect a one- parameter scaling formula for N = 
oo such as 



= t^rh 



(tlnff- 1 )= [lnGff- 1 )] 1/2 ^(tlnff- 1 ) 



(24) 

which is confirmed in Fig. 3. We numerically calculate 
to by integrating the recursion equations of g n , K n , H n , 
and their derivatives. 

Next, we investigate the property of the critical phase 
(t > 0). By considering the linear stability of Eqs. (4) 
and (5) at the points on the fixed line, we obtain 
K n+1 = e nyK ^K n and H n+1 = e nyH ^H n where K n = 
K n - K*(t), y K {t) = ln±=£, and y H {t) = ln^. 

These lead to K n oc e nyK and H n oc e nyH . While y K {t) 
is negative and therefore K is irrelevant in the critical 
phase, yn{t) is positive and approaches yd — In 2 from 
below as t ^ 0+. 

The free energy in the critical phase satisfies 

g(H,N- v ) = e~ nyd g{He nyH ^ t \N- 1 e nyd ). 

Again we omit the irrelevant parameters K and t from 
the arguments, but t-depcndcncc is included in yn- The 
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£-th order derivative of g is written as 

<7^(iJ,iV -1 ) = e n(lVH(t) - yd) g He {He nyH{t \N- 1 e nyd ). 

Since g/jj(t) < Ud, the first order derivative, i.e., the mag- 
netization, is zero for H = and N = e nyd — >• oo in the 
critical phase. The susceptibility for H — and N = e nyd 
is 

X = N'^g^ (0,1) oc JV*W (25) 

4,(t) =2 M® -1 = 1- A m (i + t ). (26) 
yd In 2 

For 1/2 < j < 1/V2, V 

varies from 1 to as shown in 
Fig. 4. In this region, \ diverges in the thermodynamic 
limit. For j > 1/V2, x is finite but higher order deriva- 
tives, i.e., nonlinear susceptibilities, diverge; the £th or- 
der derivative of the free energy with H diverges in the 
region where iyiiif) > Vd- 

As noted in Eq. (1), ip is related to the correlation 
length £, and the fact that ip approaches 1 from below 
means the divergence of £. Here we consider the finite 
size effect of this singularity. The inset of the Fig. 4 
shows the fractal exponent for size N defined as ipN = 
l°§2[x(-W)/x(-W/2)]- If we assume a power- law divergence 
of £, a finite size scaling formula, 

C 1 = 1- V>;v = t u F(t v [kiN} fl ) 

f t v for ^[lnA^ »1 

\ [IniV] M for i l/ [lniV]' 1 < 1 ' 1 J 

is expected for H = 0. If In TV is a proper cutoff length, /i 
should be unity. Equations (26) and (23) lead to 1 — 4>n oc 
t and 1 — i/jn <x [IniV] -1 , respectively. Thus v = 1 and 
/j = 1. We confirm the scaling behavior Eq. (27) in Fig. 5. 

We have investigated the 2-state Potts (Ising) model 
in the simple hierarchical network with the real-space 
RG method. The phase transition between the ordered 
and critical phases is governed by the PFBP at K = 




1 1 1 1 1 1 1 1 1 1 

0.0 0.2 0.4 . 0.6 0.8 1.0 
J 



FIG. 4: (color online) The fractal exponent is plotted with 
respect to the BBE coupling constant. The inset enlarges the 
finite size behavior i/jn around j — j c — 1/2. 




10~ 3 10~ 2 10' 10° 10' 10 2 

t (log 2 N - 0.5) 

FIG. 5: (color online) Scaling plot of the fractal exponent i/jn 
in the critical phase. Data with tpN > 0.60 are used. 



oo independently of the bare value of K. The singular 
behavior of the present model is summarized in Table I, 
where we also show the Potts model with q = l 11 and 
q > 3 15 . See also the inset of Fig. 2 for the fixed lines for 
various q systems. 

Bocttchcr et al. studied the bond percolation model, 
i.e., the 1-state Potts model, in the same graph with- 
out backbone (K = 0) 11 , and found a PF bifurcation 
of the RGFP similar to the present model (q = 2). 
They also calculate the maximum cluster size (s max ), 
which is a quantity corresponding to unconnected sus- 
ceptibility \ = X + iVm 2 because both quantities are 
defined as a summation of the two-point correlation. 
For q = 1, however, discontinuity of (s max )/iV is ob- 
served at the transition point. This is in contrast to 
the fact that x/iV oc t for q = 2. By taking the ana- 
lytic continuation of the recursion relation, correspond- 
ing to Eq. (6), for q > 3 in Ref. 15 , we obtain the PFBP 
(i>&) = (1/2,0) for q < 2 and the recursion relation at 
j = 1/2: k n+x -k n = -^k 2 -^-k 3 + 0(k 4 ). Thus the 
FP of q = 2 is marginal and has weaker stability. This is 
presumably the origin of the continuity of the transition. 

For q > 3, the saddle-node (SN) bifurcation of the FP 
is observed 15 . Consequently, two kinds of singularity ap- 
pear depending on k — e~ K : ES corresponding to SNBP 
for k > fcsN and PLS corresponding to USFP for k < A:sn- 
The SNBP (j'sn^sn) approaches the PFBP (1/2,0) as 
q 2. While the PLS and ES for q > 3 are governed by 
USFP and marginally unstable SNBP, respectively, the 
phase transition for q < 2 is governed by the stable FP. 
Thus the generalized scaling theory in Ref. 15 assuming 
instability of a FP cannot be applied to the latter. We 
emphasize that the stability or instability of the FP of 
the transition point is the most fundamental criterion of 
phase transitions. 
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q FP 
1 PFBP 



X 



N\t\' 



2 PFBP y/\t\, [XnH- 1 ]- 1 ' 2 N\t\, H' 1 



>3 SNBP 
> 3 USFP 



e i*|- 1/2 , ^ 



gsl/(l--0) 
[IniV] 2 
t _1 ,lnJV 



1/5-1 



t" 



TABLE I: Summary of singularities. We show the singular 
formula of m and x m the ordered phase t < and of £ in the 
critical phase t > 0. A discontinuous change is noted as t°. 
The fact that £ oc [IniV] 2 for q = 1 is our preliminary result 
of a Monte Carlo simulation. 
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